Universal efficiency at optimal work with Bayesian statistics 
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If the work per cycle of a quantum heat engine is averaged over an appropriate prior distribution 
for an external parameter a, the work becomes optimal at Curzon-Ahlborn efficiency. More general 
priors of the form n(a) oc 1/a'' yield optimal work at an efficiency which stays close to CA value, 
in particular near equilibrium the efficiency scales as one-half of the Carnot value. This feature is 
analogous to the one recently observed in literature for certain models of finite-time thermodynamics. 
Further, the use of Bayes' theorem implies that the work estimated with posterior probabilities also 
bears close analogy with the classical formula. These findings suggest that the notion of prior 
information can be used to reveal thermodynamic features in quantum systems, thus pointing to a 
new connection between thermodynamic behavior and the concept of information. 
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The connection between thermodynamics and the concept of information is one of the most subtle analogies in our 
physical theories. It has played a central role in the exorcism of Maxwell's demon It is also crucial to how we may 
understand and exploit quantum information • To make nanodevices [1, Q that are functional and useful, we 
need to understand their performance with regard to heat dissipation and optimal information processing. To model 
such systems, standard thermodynamic processes and heat cycles have been generalised using quantum systems as the 
working media It is well accepted that the maximal efficiency, rye = 1 — T2/T1, where Ti(T2) is the hot (cold) 

bath temperature, is only obtained by a reversible heat engine which however involves infinitely slow processes. For 
heat cycles running in a finite time, the concept of power output becomes meaningful. Curzon and Ahlborn JTTj , first 
of all displayed an elegant formula in the so called endoreversible approximation for efficiency at maximum power, 
yycA = 1 — ~ Vc- The appearance of an optimal efficiency in different models with a value close to Curzon-Ahlborn 
(CA) value, has raised the issue of its universality that has captured the imagination of workers in this area since many 
years [H]. More recently [il],[i3-[23, a universal form for optimal efficiency, has been discussed within finite-time 



thermodynamics in the near-equilibrium regime (small value of 77c) as given by 77 ~ '7c/2 + 0(rjc^). 

On the other hand, in recent years, Bayesian methods of statistical inference have gained popularity in physics [23 |. 
In Bayesian probability theory, the central role is played by the concept of prior information. It represents state of 
our knowledge about a system before any experimental data is acquired. The assignment of a unique distribution to 
a given prior information is a non-trivial issue but may be argued on the basis of maximum entropy principle and 
certain requirements of invariance [26l. [27j. Using Bayes' theorem [25| one can then update the prior probabilities 
based on the new information gathered from the data. Recently, Bayesian methods have been applied to modular 
structure of networks [2^ , inference of density of states [29i] , the interpretation of quantum probabilities [s^ and other 
inverse problems [sij . 

In this letter, a Bayesian approach is used to show that the efficiency at optimal work for a quantum heat engine 
is related to CA value, after the work per cycle is averaged over the prior distribution of an external parameter. In 
contrast to the finite-time models (23] , the heat cycle considered here is performed infinitely slowly. The application 
of Bayes' theorem gives the optimal efficiency exactly at CA value for a whole class of priors and for arbitrary bath 
temperatures. The present analysis thus provides a novel argument for the emergence of thermodynamic behavior in 
quantum heat engines from a Bayesian perspective. 

As a model of a heat engine, consider a quantum system with Hamiltonian Hi = X^n I") (^L where energy 

eigenvalues Sn ^ = £nai. The factor e„ depends on the energy level n as well as other fixed parameters/constants 
of the system; ai is a controllable external parameter equivalent to say, the applied magnetic field for a spin-1/2 
system. Other examples of this class are I-d quantum harmonic oscillator (ai equivalent to frequency) and a particle 
in 1-d box (ai inversely proportional to the square of box- width). Initially, the quantum system is in thermal state 
p{ai) = Pn^^ 1'^) ('^1 s-t temperature Ti with its eigenvalues given by the canonical probabilities Pn^ . The quantum 
analogue of a classical Otto cycle between two heat baths at temperatures Ti and T2 involves the following steps 
[1^ : (i) the system is detached from the hot bath and made to undergo the first quantum adiabatic process, during 

('2') (2) 

which the system hamiltonian changes to H2 = where e„ = £1102, without any transitions between the 
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levels and so the system continues to occupy its initial state. For 02 < oi, this process is the analogue of an adiabatic 
expansion. The work done 6j/ the system in this stage is defined as the change in mean energy Wi — TT(p(ai)[H2—Hi])] 

(2) 

(ii) the system with modified energy spectrum e„ is brought in thermal contact with the cold bath and it achieves 

a thermal state p{a2) — X^nP" The modified canonical probabilities pii now correspond to temperature T2. 

On average, heat rejected to the bath in this stage is defined as Q2 — Tr([p(a2) — p{ai)]H2); (iii) the system is now 
detached from the cold bath and made to undergo a second quantum adiabatic process (compression) during which 
the hamiltonian changes back to Hi. Work done on the system in this stage is >V2 = Tr(p(a2)[-ffi — -^2]); (iv) finally, 
the system is brought in thermal contact with the hot bath again. Heat is absorbed by the system in this stage whence 
it recovers its initial state and its temperature attains back the value Ti. The total work done on average in a cycle 
is calculated to be 



n 

= (ai-a2)E£„(p«-p(2)) >0. 



(1) 
(2) 



Similarly, heat exchanged with hot bath in stage (iv) is given by Qi — aiJ2n^n (^PrP — Pr?^ > 0. Heat exchanged 
by the system with the cold bath is Q2 = W — Qi < 0- The efficiency of the engine 77 = W/Qi, is given by 



1 I "2 

7? = 1 - — 



(3) 



For convenience, we express W = >V(ai, ry), using Eq. Consider an ensemble of such systems where now the value 
of parameter oi may vary from system to system. If the ensemble corresponds to an actual preparation according to 
a certain probability distribution H(ai), then the state of the system can be expressed as p = J p{ai)Il{ai)dai. Each 
system in the ensemble is made to perform the quantum heat cycle described above, with a fixed efficiency 77. We 
wish to study the optimal characteristics of the average work, in particular the efficiency at which the work becomes 
optimal. Clearly, choice of the probability distribution H(ai) is expected to play a significant role in the conclusions. 
In the following, we analyse this problem by choosing a distribution H(ai) according to the prior information available 
and show that the efficiency at optimal work is closely associated with CA value. 

For simplicity, we now consider a two- level system as our working medium, with Eq = and ei = 1, so that the 
initial energy levels are and ai . The work over a cycle in this case is 



W(ai,?7) = aiT] 



1 



1 



(l + e''i/^i) (l + e°i(i-")/^2) 



>0, 



(4) 



where Boltzmann's constant is put equal to unity. The average work with the initial state p for a given 77, can be 
expressed as 



W 



W{ai,ri)Il{ai)dai. 



(5) 



A central issue in Bayesian probability is to assign a unique prior distribution corresponding to a given prior informa- 
tion. If the only prior information about the continuous parameter oi is that it takes positive real values but otherwise 
we have complete ignorance about it, then Jeffreys has suggested the prior distribution H(ai) oc l/oi [H, |3l|, or in 
a finite range, H(ai) = [In (amax/amin)] ^ (l/ai)j where Omin and Cmax are the minimal and the maximal energy 
splitting achievable for the two-level system. For the above choice, we obtain 



W = 



In 



1 + ea„,x(l-»?)/T2 

(1 - r]) V 1 + ea„.n(i-»;)/T. 



■In 



-Tiln 



1 



1 + 



(6) 



It can be seen that the average work W vanishes for i] — and i] = rjc- In between these values of j], the average work 
exhibits a maximum. We look for the efficiency at which this work becomes maximal for the given range [flmin, amax], 
by imposing the condition dW/djj = 0. Here we consider the limit of a^nin which gives 



Tn 



{i-r,y 



■In 



I _,_ gQ„,ax(l-»))/T2 



Tiln 



1 



V 



= 0. 



(7) 
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The solution 77 of this equation has been plotted against 
flmax » Ti, the above expression reduces to 



Ti 



To 



in Fig. 1. Interestingly, in the asymptotic limit of 



(8) 



which yields the efficiency at optimal work as 77 = 1 — a/ T2 / Ti , exactly the C A value. More significantly, the conclusion 

also holds in general i.e. for a working system with spectrum Sn'^ = and with Jeffreys' prior. It is to be noted 
that in the asymptotic limits, the expression for average work (Eq. ([6])) diverges. However, the limits are taken after 
the derivative of work is set equal to zero in order to obtain well-defined expressions for the efficiency. 




FIG. 1: Efficiency versus flmax using Eq. 0. The curves correspond to T2 = 1 and Ti taking values 2,4,6 respectvely, from 
bottom to top. Apart from the approach to corresponding CA value at large Omax, it is also seen that the limit is approached 
slowly for larger temperature differences. 



It is conceivable that other choices of the prior may yield similar results. To study consequences of deviations from 
the above choice, we consider a class of prior distributions, n(ai) = iVai~^, defined in the range [0, Omaxlj where 
= (1 — 7)/amax^~''' and 7 < I. Upon optimisation of the average work as defined in Eq. ([5]) over 77, we get 



(ai) 



1-7 



(ai) 



1-7 



l + e°i/'^i 1 + eai(i-r))/T2 



da^ — -j- 

T2 



(a^)2-7eai(l-r^ 

(1 + e°i(i-'?)/^2)2 



dai — 0. 



(9) 



In the limit a,nax becoming very large, the above integrals can be evaluated using the standard results ^] . Then the 
above equation is simplified to 



32-7 



0, 



(10) 



where ~ T2/T1. Now as 7 ^ 1, the above equation reduces to Eq. ([5]) and so CA value is also a limiting value for 
this model. Interestingly, even for other allowed values of 7, the solution rj* of Eq. (|10p depends only on the ratio 9, 
apart from the parameter 7. In particular, Laplace and Bayes have advocated a uniform prior to quantify the state 

(1 + ^*)02 ^ 0, 



of complete ignorance. For this case, we set 7 = 0. Then the above equation becomes (1 — rj 

1/3 



which has only one real solution given as 

77* = 1 + 



94/3 



3 1 



1 -I- — 



1/3 




(11) 



This solution along with other numerical solutions of ([TU]) for general 7 < 1 are shown in Fig. 2. Remarkably, these 
curves stay very close to the CA value. However, at this point it is not possible to say in general what prior information 
may be quantified by the parameter 7. The curves in Fig. 2 are also closely similar to those observed in finite-time 
models at optimal power 2^. It is seen here that in the near-equilibrium regime, all the curves merge into each other 
and approach the CA value which is approximately 77^/2 in this limit. This can be shown as follows: taking 6 to be 
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close to unity in the near-equilibrium case, rjc = {1 — 6) is close to zero. The efficiency r]* being bounded from above 
by the Carnot value is thus small too. On using these facts in the expansion of Eq. ([TU]) . we get 



V 



Vc 



(3-7 ). 
16 



(12) 



Thus we recover the linear term ?7c/2 mentioned earlier. For general values of 9, the CA value is a lower bound for 
the efficiency at optimal work when < 7 < 1 [13] ■ 
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FIG. 2: Optimal efficiency rj* versus 6 — T2/T1, with parameter 7 of the prior distribution taking values 0,0.25,0.75,1.50 
respectively, from top to bottom, excluding the solid line which represents the CA values, 1 — y/6. For close to equilibrium (9 
nearly unity), the optimal efficiency exhibits a universal form, independent of 7 and given by ?;* « (1 — 9)/2 — rjc/2. 



So far we have observed that the use of Jeffreys' prior implies that efficiency at optimal work approaches CA value 
for arbitrary bath temperatures. Further, the efficiency also approaches a universal form for a class of priors, for nearly 
equal bath temperatures. In the following, we show that application of Bayes' theorem can restore the efficiency back 
to the exact CA value even for the latter choice. Bayes' theorem gives a prescription to convert the prior probabilities 
n(ai)dai into posterior probabilities. Note that during the ffist quantum adiabatic process on the two-level system, 
the energy levels change from (0, ai) to (0,02), but the system continues to occupy its initial state. The respective 
occupation probabilities are now interpreted as conditional probabilities, given by p(t \ai) = 1/(1 + exp(ai/Ti)) and 
p{], \ai) — l/(H-exp(— ai/Ti)). If the system is found in the up (t) state, the work done in this step is (02 — 01) = — oi?? 
and the posterior probabilities are given by 



p{ai \ t)dai = 



pit |ai)n(ai)dai 
/ pit \ai)U{ai)dai' 



(13) 



The average work for this process is now given by Wi — J {—air])p{ai\ t)rfai- On the other hand, if the system is 
found in the down {\.) state, the work is zero. Similarly, for the second quantum adiabatic process, the work performed 
can be either {+airi) or and the avarage work W2 for that process can be similarly calculated using the respective 
posterior probabilities. Now choosing the prior n(ai) — Nai^"' and with the system being in up state, the average 
work for the total cycle [Wi + W2), in the limit of large Cmax is given by 



W(?7) 



(2'^-i-l)(l~7)C[2-7] 
(27-l)C[l-7] 



(1-^) 



(14) 



where 7 < 1. So using posterior probabilities, a well defined expression for average work is obtained even if the prior 
is non-normalisable in the asymptotic limit. More generally, given the value of external parameter ai and assuming 
canonical probabilities p{n\ai) to find the system in nth state, we infer the probability p{ai\n)dai about the value of 
ai, if the system is actually found in the nth state. Remarkably, the work given by eq. (ITU) attains optimal value 
exactly at the CA efficiency, regardless of the value of 7 in the prior. Furthermore, the average work W{ri) shows the 
same dependence on efficiency 77 as found for the classical Otto cycle in fisf . 

In conclusion, we have argued the emergence of CA value as the efficiency at optimal work in quantum heat engines 
within a Bayesian framework. This effect of incorporating Bayesian probabilities leading to classical thermodynamic 
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behavior in quantum systems has not been addressed before and may shed new hght on the connection between 
information and thermodynamics. Due to current interest in smah scale engines, the observation of similar curves 
(Fig. 2) as obtained in some recently proposed models of these engines, points to an interesting link between finite 
time models and our model based on the idea of prior information. Addressing these issues would hopefully lead to 
a broader perspective on the performance characteristics of small engines and also help to understand the limits of 
their performance based on principles of information. 
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